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ABSTRACT 

<N ■ 

^ ' Context. 3C 279, the first quasar discovered to emit VHE y-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 
CNj j during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high 
. y-ray state. 

l/~) • Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2<r, corresponding to a F(> 150 GeV) (3.8 ± 0.8) • 
C^l ' 10~ u ph cirT 2 s while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the 
* \ source. We study the multiwavelength behaviour of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical 
, and X-ray frequencies and for 2009 also y-ray data from Fermi. 
I*— *j ■ Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs 
j ' (including the February 2006, which has been previously published) are modelled with one-zone inverse Compton models and the emission on 
. . ' January 16, 2007 also with two zone model and with a lepto-hadronic model. 

^ . Results. We find that the VHE y-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a 
• ^ ■ jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily. 

% , ■ Key words, gamma rays: galaxies — quasars: individual: 3C 279 



1. Introduction 

3C 279 was the first quasar discovered as a y-ray source with the 
Compton Gamma-Ray Observatory (Hartman et al. 1992), and 
the first flat-spectrum radio quasar (FSRQ) discovered to emit 
very high energy (VHE, defined here > 100 GeV) y-rays (Albert 
et al. 2008a). With a redshift of 0.536 (Hewitt & Burbridge 
1993), 3C 279 is also the most distant of the VHE y-ray emitting 
sources discovered so far. 

3C 279 is one of the brightest quasars at all wavelengths and 
its multiwavelength behavior and jet structure has been studied 
in great detail (e.g. Hartman et al. 2001, Chatterjee et al. 2008). 
Its relativistic jet, which is the source of the radio to VHE y- 



ray emission, is closely aligned with the line of sight (the angle 
varies, but is sometimes as small as < 0.5°, Jorstad et al. 2004). 
The radio to optical emission is synchrotron radiation emitted 
by the relativistic electrons spiraling in the magnetic field of 
the jet. In this low-energy regime the total flux density varia- 
tions are well described by shocks propagating in the jet (e.g. 
Lindfors et al. 2006). The X-ray emission can be explained by 
the synchrotron self-Compton mechanism (SSC, e.g. Maraschi 
et al. 1992, Hartman et al. 2001, Sikora et al. 2001), where the 
synchrotron photons emitted by the jet act as seed photons for 
the inverse Compton scattering. 
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However, there is no consensus about the emission mecha- 
nism and site of the y-ray and VHE y-ray emission in 3C 279. 
The emission can, in principle, be explained by both leptonic 
and hadronic models: the leptonic models mostly rely on ex- 
ternal Compton (EC, e.g. Hartman et al. 2001, Albert et al. 
2008a), invoking the inverse Compton scattering of external 
photons from accretion disk or broad line region (BLR) clouds 
(Dermer & Schlickeiser 1993, Sikora et al. 1994) , while in the 
hadronic models (Mannheim & Biermann 1992, Miicke et al. 
2003, Bottcher et al. 2009) the VHE y-ray photons are pro- 
duced by proton initiated cascades or directly through proton 
synchrotron radiation (but see Sikora et al. 2009 for criticisms). 

The leptonic models are very sensitive to the site of the emis- 
sion: the external Compton models relying on photons originat- 
ing from broad line emission clouds are not efficient if the emit- 
ting blob is outside the BLR and in the SSC models the y-ray 
emission must originate from a different emission region than 
the main component of the synchrotron radiation in order to re- 
produce the observed y-ray flux (Bottcher et al. 2009). It should 
also be noted, that independently of the emission mechanism the 
internal absorption cannot be neglected if the emission region 
is located inside the BLR (Sitarek & Bednarek 2008, Tavecchio 
& Mazin 2009). Alternatively, the emission can be produced in 
regions located far beyond the broad line region, at distances at 
which the dominant radiation field for EC is that of the parsec- 
scale dusty torus (Sikora et al. 2008). In this case, internal ab- 
sorption can be neglected up to ~ 1 TeV but, due to the large 
size of the emission region, a minimum variability timescale of 
the order of ~ 1 day is expected. 

Amongst Fermi detected blazars (Abdo et al. 2009) hard 
overall spectrum (spectral index=2.32) but also the compara- 
tively weak evidence of a break in the spectrum at 1 GeV. The 
spectral index after the break energy is given as 2.50 without ap- 
parent cut-off. This is one of the hardest spectrum of all bright 
FSRQs above a few GeV and makes it one of the prime target 
for VHE y-ray observations. 

In this paper we present VHE y-ray observations of 3C 279 
performed by the MAGIC-I telescope in January 2007 and from 
December 2008 to April 2009. The 2007 observations were trig- 
gered by an optical outburst in the source while the December 
2008 observations were triggered by an alert from the Fermi 
space telescope. The January-April 2009 data were taken as a 
part of the followed multiwavelength campaign. We present the 
data analysis and results of the MAGIC observations, together 
with simultaneous multiwavelength observations and a discus- 
sion of theoretical models. For comparison and completeness 
we also summarize the 2006 observation campaign (Albert et 
al. 2008a). The plan of the paper is as follows. The MAGIC 
observations, data analysis and results are presented in Section 
2. Multiwavelength observations are presented in Section 3 and 
discussed in Section 4 together with the MAGIC data. VHE 
y-ray and multifrequency data are combined in Section 5 to 
build quasi-simultaneous spectral energy distributions at differ- 
ent epochs and theoretical models are presented and critically 
discussed. A summary of results and conclusions are given in 
Section 6. 

2. MAGIC Observations 

The data described in this paper were taken with the first MAGIC 
(Major Atmospheric Gamma-ray Imaging Cherenkov) telescope 
as a standalone imaging atmospheric Cherenkov telescope. It is 
located on the Canary Island of La Palma. MAGIC has a stan- 
dard trigger threshold of 60 GeV for observations at low zenith 



angles, an angular resolution of ~ 0.1° on the event by event 
basis and an energy resolution above 150 GeV of ~ 25% (see 
Albert et al. 2008b for details). 

MAGIC observed 3C 279 in January 2007 and from 
December 2008 to April 2009. Due to changes in the telescope 
systems these data sets were analyzed separately. 

2.1. January 2007 

3C 279 was observed during nine nights in January 2007 for a 
total of 23.6 h, 18.6 h (seven nights) of which passed the quality 
selection. The observations were conducted in moon, dark night 
and twilight. The zenith angle range of the observations was 35 
to 46 degrees. The readout chain included a fast 300 MSample/s 
FADC system, which allows us to use the time information of 
the showers in the analysis as described in Aliu et al. (2009). 
Since the average night sky background noise level stayed below 
3 times the dark time level, all data has been analyzed with the 
same standard image cleaning levels as discussed in Britzger et 
al. (2009). The data were taken in On mode where the telescope 
is pointing directly towards the source. So-called Off data were 
taken under similar conditions (zenith angle, trigger rates and 
weather conditions) from September 2006 until January 2007 
and has been used to estimate the background level. Due to the 
installation of the new 2 GHz MUX-FADC read-out system in 
2007 February, no Off data after January 2007 were fulfilling 
the requirement of similar observation and hardware conditions. 

The 2007 data were analyzed using so-called Wiirzburg 
Analysis chain described in detail in Bretz & Wagner (2003) and 
Bretz & Dorner (2008). The data were calibrated following the 
description in Albert et al. (2008c), the signal was extracted at 
the pulse maximum using a spline method, the air-shower im- 
ages were cleaned of noise from night-sky background light by 
applying a time image cleaning. In the first step a minimum num- 
ber of 6 photo-electrons in two adjacent pixels is required (so 
called "core pixels"). Each pixel next to a core pixel, which is 
above a threshold of 3 photo-electrons is considered a "boundary 
pixel". Additionally to these limits in the charge of the pixels it 
is required that the arrival time of core pixels is within 1 .75 ns of 
the mean arrival time of all core pixels and that boundary pixels 
arrive within 4.5 ns of the mean arrival time of the core pixels. 
Both the charge and the time limits have to be fulfilled, other- 
wise the pixel is considered to contain noise and deleted from 
the image. We parametrized the shower images (Hillas, 1985) 
and used a SIZE-dependent parabolic cut in AREA, WIDTH 
and LENGTH (Riegel et al. 2005) for the y/hadron separation. 
The cut parameters were optimized on Crab Nebula data taken 
within the same zenith angle range as the 3C 279 data. The en- 
ergy threshold of the analysis was 220 GeV. The arrival direction 
of the y-rays is reconstructed using the DISP method (Fomin et 
al., 1994; Lessard et al., 2001), which was adapted to use the 
shower timing information as described in Aleksic et al. (2010). 
The significance of a detection is evaluated by comparison of 
the number of events in the On and the Off data sample using 
formula 17 of Li & Ma (1983). 

Since in 2006 the source was detected in a single day flare, 
every observation day was analyzed separately. Table 1 summa- 
rizes the analysis results of the January 2007 campaign. Out of 
the seven observation nights with good data quality only one, 
MJD 54116 (January 16), revealed a significant excess in the 
MAGIC data after Off-subtraction. The significance of the ex- 
cess is 5.6cr (pre-trial) with 64.1 excess events above 102.9 nor- 
malized Off events in 151 min of observations. The 9 2 (squared 
distance between the true and reconstructed source position, see 



2 



Aleksic et al.: 3C 279 

Table 1. Night-by-night results of MAGIC January 2007 observations of 3C 279. 



Observation night [MJD] 


Observation time [min] 


Excess events [counts] 


background events [counts] 


On/Off scaling factor 1 


significance 2 


54115 


149.5 


17.6 


97.4 


0.058 


1.7 


54116 


151.1 


64.1 


102.9 


0.062 


5.6 


54117 


157.2 


-7.9 


100.9 


0.060 


-0.8 


54118 


153.6 


-6.0 


82.0 


0.049 


-0.7 


54120 


164.1 


-12.3 


117.3 


0.07 


-1.1 


54121 


166.7 


12.7 


119.3 


0.07 


1.1 


54122 


175.1 


18.2 


126.8 


0.076 


1.5 



1 The scaling factor is the ratio between the on and the off distribution normalized outside the signal region. 

Significance is given as standard deviations and is not corrected for seven trials, arising from diving the data sample by seven 
(corresponding to seven nights with good quality data). 
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Fig. 1. Reconstructed shower direction 2 for the Off (gray 
shaded area) and On (black crosses) data as observed on January 
16, 2007 by MAGIC. The vertical dotted line corresponds to the 
apriori defined signal region. An excess of 64 events is clearly 
visible in the On data as compared to the normalized Off data. 
The corresponding pre-trial significance is 5.6cr and 5.4<x after 
correction for seven trials, respectively. 



e.g. Daum et al. 1997) distribution is shown in Fig. 1. None of 
the other observation nights yielded in a significant excess. 

For spectrum (Figure 2) and light curve softer cuts that have 
a higher y-ray efficiency were used. The VHE y-ray spectrum 
of the flare can be described by a simple power law (with the 
differential flux given in units of TeV~' irT 2 s ): 



dF 



— = (5.7 + 1.3) x 10~ 7 ( ) 



•3.1±1.1 



300 GeV ^ 



(1) 



The highest spectral energy point has a medium energy of 
350 GeV. For a Crab like spectrum we have a systematic en- 
ergy scale error of 16%, a systematic error of 11% on the flux 
normalization (without the energy scale error), and a systematic 
slope error of 0.2 (a detailed list of all the contributions can be 
found in Albert et al. 2008b), while for the soft spectrum like 
3C 279 the systematic errors are expected to be slightly larger. 
Assuming that 3C 279 always emits y-ray s above 150 GeV we 
calculated the corresponding light curve (see Section 4) on the 
night by night basis. For January 16th the integral flux above 
150 GeV is (3.8 + 0.8) ■ 10" 10 ph ctrT V 1 . 

The results were cross-checked using the standard MAGIC 
analysis chain (Albert et al. 2008c). The standard analysis re- 
sulted in lower significance, 4.6cr, for the detection. For flux 
and spectral index the results were comparable for both analy- 
sis chains. 
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Fig. 2. Spectrum of the VHE y-ray flare January 16, 2007. The 
blue points correspond to the measured spectrum with the re- 
spective errors on the flux and the bin size in energy. The re- 
sults of a simple power law fit are shown as blue line. The gray 
shaded region represents the systematic error of the analysis. The 
red points show the spectrum deabsorbed using the EBL model 
by Dominguez et al. (2010) and the red dashed line the fit to 
this EBL corrected spectrum. For comparison the spectrum from 
2006 (Albert et al. 2008a) is shown with green points. 



2.2. December 2008- April 2009 

MAGIC observations started on December 9, 2008 after the 
Fermi Collaboration announced a high GeV y-ray state of the 
source (Ciprini et al. 2008). Unfortunately 3C 279 could only 
be observed at zenith angles larger than 46 degrees and for a 
short time (~ 10 min) due to visibility constraints. These short 
runs were insufficient to detect the source. Follow-up observa- 
tions were conducted at the end of December 2008 until April 
2009 under more favorable zenith angles, smaller than 35 de- 
grees. 

A total observation time of 28.1 h was accumulated over 20 
days with the 2 GHz MUX-FADC read-out system. 35 min were 
recorded in On-mode, while the remaining data (27.5 h) were 
taken in wobble mode where the source was displaced by 0.4 ° 
from the camera center in order to allow the recording of si- 
multaneous Off data with the same offset from the camera cen- 
ter (Daum et al. 1997). After quality selection 1 1 .9 h of data 
were used in the analysis. The main part of this data set is from 
January 2009, but it also includes one night (29) in December 
2008 and one (16) in April 2009. 

The data were analyzed using the standard analysis chain 
as described in Albert et al. (2008b, 2008c) and Aliu et al. 
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Table 2. Night-by-night results of MAGIC December 2008 - April 2009 observations. 



Observation night [MJD] 


Observation time [min] 


Excess events [counts] 


background events [counts] 


On/Off scaling factor 1 


significance' 1 


54829 


50.4 


-0.3 


65.3 


0.03 


-0.04 


54832 


60.0 


-10.7 


92.7 


0.02 


-0.98 


54834 


81.0 


4.8 


103.2 


0.02 


0.41 


54838 


78.0 


-1.4 


48.4 


0.08 


-0.18 


54852 


33.6 


12.9 


42.1 


0.07 


1.73 


54855 


157.8 


4.1 


241.9 


0.06 


0.24 


54856 


140.4 


5.3 


70.7 


0.10 


0.57 


54937 


114.6 


-0.3 


110.3 


0.03 


-0.03 



1 The scaling factor is the ratio between the on and the off 6 distribution normalized outside the signal region. 

Significance is given as standard deviations and is not corrected for trials, arising from diving the data sample by eight 

(corresponding to eight nights with good quality data). 



(2009). In order to suppress the background showers produced 
by charged cosmic rays, a multivariate classification method 
known as Random Forest is used (Albert et al. 2008d). For 
every event, the algorithm takes as input a set of image pa- 
rameters, and produces one single parameter as output, called 
Hadronness. The background rejection is then achieved by a cut 
in Hadronness, which was optimized using Crab Nebula data 
taken under comparable conditions. 

A cut in 9 2 , was used to extract the signal and also optimized 
in the same way. An additional cut removed the events with a 
total charge of less than 200 photo-electrons (phe) in order to 
assure a better background rejection. 

We find no signal in this data set. The number of excess 
events is 29 over a background of 775 events, corresponding to 
a significance of 0.94<x. The night by night results are reported 
in Table 2. 



Table 3. Differential upper limits (95% confidence level) for 
December 2008-April 2009 observations. 



Energy bin[GeV] 


Upper Limit[l/TeV/m2/s] 


150-200 


1.810 3 


200-315 


1.6-10- 6 


315-430 


7.2-10-' 


430-580 


2.5-10-' 



term optical R-band light curve covering all the MAGIC obser- 
vation periods discussed in this paper. In 2009 also polarimetric 
observations of 3C 279 were performed. The photometric optical 
R-band observations are performed with the Tuorla 1 meter tele- 
scope (Finland) and Kungliga Vetenskapsakademien (KVA) 35 
cm telescope (La Palma). The latter can be controlled remotely 
from Tuorla Observatory. 
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Fig. 3. Long term optical R-band observations of the Tuorla 
blazar monitoring program of 3C 279. The giant flare near the 
beginning of 2007 is clearly visible. 



All data were combined for the calculation of flux upper lim- 
its in the energy range between 150 GeV and 580 GeV (higher 
energies can be neglected due to the strong EBL absorption). 
Including a 30% systematic error and assuming a power law 
with a spectral index of -4 the integrated flux upper limit (2cr 
c.l.) is 6.3T0 -7 ph/m 2 /TeV/s. In Table 3 the differential values 
are reported. 

3. Multiwavelength Observations 

We present here long term optical monitoring data from the 
Tuorla blazar monitoring program, and near infrared data from 
REM (for January 2007), X-ray data from SWIFT (for February 
1, 2009) and y-ray data from Fermi (for January 2009). 

3.1. Optical Observations 

3C 279 has been observed regularly since 2004 as a part of the 
Tuorla blazar monitoring program^ Figure 3 shows the long 

1 http://users.utu.fi/kani/lm/ 



KVA consists of two telescopes, the larger one being a 60cm 
(f/15) Cassegrain telescope equipped with a CCD polarimeter 
capable of polarimetric measurements in BVRI-bands using a 
plane-parallel calcite plate and a super-achromatic lambda/2 re- 
tarded For 3C 279 the observations were done without filter to 
improve the signal-to-noise of the observations. 

The 35cm telescope is used for photometric measurements in 
B, V and R-bands. The photometric measurements are made in 
differential mode, i.e. by obtaining CCD images of the target and 
calibrated comparison stars in the same field of view (Raiteri et 
al. 1998). The magnitudes of the source and comparison stars are 
measured using aperture photometry and the (color corrected) 
zero point of the image determined from the comparison star 
magnitudes. Finally, the object magnitude is computed using the 
zero point and a filter-dependent color correction. Magnitudes 
are then transferred to linear fluxes using the formula F = Fq ■ 
10 mag ^ 2 5 , where mag is the object magnitude and is a filter- 
dependent zero point (in the R-band the value Fy = 3080 Jy is 
used from Bessell et al. 1979). 

For polarimetric measurements the normalized Stokes pa- 
rameters and the degree of polarization and position angle were 
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calculated from the intensity ratios of the ordinary and extraor- 
dinary beams using standard formula (e.g. Degl'Innocenti et al. 
2007) and semiautomatic software specially developed for polar- 
ization monitoring purposes. During some of the nights, polar- 
ized standard stars from Turnshek et al. (1990) were observed to 
determine the zero point of the position angle. The instrumental 
polarization of the telescope has been found to be negligible. 

3.2. REM observations in January 2007 

In January 2007 3C 279 was intensively monitored by REM 
in the IR band (J, H and K filters). REM is a 60-cm diameter 
fast-reacting telescope with 10° s _1 pointing speed located at the 
Cerro La Silla premises of the European Southern Observatory, 
Chile (Zerbi et al. 2001; Chincarini et al. 2003; Covino et al. 
2004a; Covino et al. 2004b). The telescope hosts REMIR, an 
infrared imaging camera, and ROSS, an optical imager and slit- 
less spectrograph. The two cameras observe simultaneously the 
same field of view of 10'xlO' thanks to a dichroic mirror, al- 
though only the infrared data are reported here. The telescope 
was designed to rapidly follow transient high-energy events as 
Gamma-Ray Bursts (GRBs) but in the idle time it is perfectly 
suited for monitoring programs of variable sources. The obser- 
vations, including the short and long-term scheduling, are totally 
automatic. Typical exposure durations were of 150 s in the /, H 
and Ks filters. Data were reduced in a standard way by means 
of tools provided by the ESO-Eclipse package (Devillard 1997). 
Standard aperture photometry was derived and results calibrated 
by a suitable number of well-exposed 2MASS objects in the 
fieldQ 

3.3. Swift and Fermi observations in January-February 2009 

The Swift satellite (Gehrels et al. 2004) is a NASA mission, 
launched in 2003. Swift performed Target of Opportunity obser- 
vation of 3C 279 on February 1 , 2009. 

Swift/XRT is a Wolter type I grazing incidence telescope, 
with 110 cm 2 effective area, 23.6' arcminutes FOV and 15" an- 
gular resolution, sensitive in the 0.2-10 keV energy band. Data 
were reduced using the software distributed with the heasoft 
6.3.2 package by the NASA High Energy Astrophysics Archive 
Research Center (HEASARC). The xrtpipeline was set for 
the photon counting or window timing modes and having se- 
lected single pixel events (grade 0). Data shown in the SEDs 
were rebinned in order to have at least 30 counts per energy bin. 
Power law models have been fitted to the spectra. The X-ray red- 
dening due to absorbing systems along the light travel path has 
been corrected assuming the Galactic value for column density 
of neutral hydrogen N H = 2.1 x 10 20 cm" 2 (Kalberla et al. 2005). 

UVOT is a 30 cm diffraction limited optical-UV telescope, 
equipped with 6 different filters, sensitive in the 1700-6500 A 
wavelenght range, in a 17' x 17' FOV. During the pointing of 
February 1, 2009 used below for the SED, only filters UVW1 
and U were available. Analysis was performed by means of the 
uvotimsum and uvotsource tasks with a source region of 5", 
while the background was extracted from a source-free circular 
region with radius equal to 50" (it was not possible to use an 
annular region, because of a nearby source). The extracted vF v 
magnitudes have been corrected for Galactic extinction using the 
values of Schiegel et al. (1998) and applying the formulae by 
Pei (1992) for the UV filters, and finally converted into fluxes 
following Poole et al. (2008). 

2 http://www.ipac.caltech.edu/2mass/ 



The Fermi LAT is a pair conversion telescope designed to 
cover the energy band from 20 MeV to greater than 300 GeV 
which operates in all-sky scanning mode. 

The average LAT spectrum was derived in the period of 
the MAGIC pointings using the publicly available LAT dataQ 
The photons of class 3 (DIFFUSE) with energy in the range 
0.1-100 GeV collected from January 21, 2009 to January 
31, 2009 were selected. These data were processed by using 
Science Tools 9.15.2, which includes the Galactic diffuse and 
isotropic background and the Instrument Response Function IRF 
P6_V3_DIFFUSE. Then photons in the good-time intervals and 
within a region of interest (ROI) with radius of 10° from the 
source radio position were selected and a cut on the zenith angle 
parameter (< 105°) to avoid the Earth albedo was applied. The 
following steps were to calculate the live-time cube, the expo- 
sure map and the diffuse response. 

With all these information at hands, an analysis by using an 
unbinned likelihood algorithm (gtlike) in separate energy bins 
was performed. The model included the isotropic and Galactic 
diffuse backgrounds, the source of interest, all the 1FGL sources 
in the ROI and, possibly, additional sources not included there 
but identified in the map. For all the point sources we assumed 
a power law spectrum, with flux and photon index as a free pa- 
rameter and calculated the corresponding test statistic (TS , see 
Mattox et al. (1996) for a definition; in practice one assumes 
y/TS - cr, the significance of the detection). In the SED we re- 
port the four bins with high significance, TS > 25. 

4. Multiwavelength behavior 

In addition to multiwavelength data discussed in the previous 
section, there are also previously published multiwavelength 
data for all the MAGIC observing epochs (spring 2006: Bottcher 
et al. 2007, Collmar et al. 2010, winter 2006-2007: Larionov et 
al. 2008, 1996-2007: Chatterjee et al. 2008, winter 2008-spring 
2009: Abdo et al. 2010). In the following we discuss the multi- 
wavelength behaviour of the source in these three epochs com- 
paring the data described in this paper to previously published 
results. 

In the spring of 2006 when 3C 279 was first discovered 
by MAGIC, the source was in a high optical and X-ray state. 
However no optical or X-ray flare simultaneous to the VHE y- 
ray flare was observed (Fig. 4). We are not aware of any optical 
polarization monitoring data for this epoch. In the radio bands 
the source was in quiescent state (Bottcher et al. 2007, Collmar et 
al. 2010) and there was no coincident very long baseline interfer- 
ometry (VLBI) knot emerging from the VLBI core (Chatterjee 
et el. 2008), suggesting that the VHE y-ray emission originates 
from region that is opaque to radio frequencies. 

In January 2007 3C 279 was very bright in the optical band, 
reaching a peak flux of 19 mJy (R~13) on MJD 54114 . This 
is the maximum flux observed from the source in five years of 
Tuorla blazar monitoring. However, during the MAGIC obser- 
vation the optical flux was already decreasing and on the night 
of January 16, 2007 (MJD 54116) no increase in the optical 
flux was detected. A similar behavior was seen in the near- 
IR band, with a clear peak at January 13 followed by a de- 
cay. Unfortunately REM did not observe 3C 279 on January 15 
and 16, but both on January 14 and 17 the data show a steep 
spectrum, smoothly joining with the R-band measure of KVA. 
Therefore, it is reasonable to assume that also the IR flux de- 
creased smoothly during the VHE y-ray flare. 

3 accessible from http://fermi.gsfc.nasa.gov 
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Fig. 4. Light curve in different energy bands during the 2006 
observation campaign. From top to bottom: VHE y-ray flux 
between 100-500 GeV as measured by MAGIC (Albert et al. 
2008a), RXTE PCA flux in the 2-10 keV range (Chatterjee et al. 
2008) and KVA R-band observations. 



The fluxes inferred from the MAGIC (> 150 GeV), the RXTE 
PCA (2-10 keV), the KVA (R-band) and REM near infrared ob- 
servations of January 2007 are compared as a function of time 
(light curves) in Fig. 5. As can be seen in Fig. 5, the VHE y-ray 
flare does not coincide with neither the maximum of the opti- 
cal flux or the maximum of the X-ray flux that were both ob- 
served prior to MAGIC observations. A comparison of the dif- 
ferent light curves suggests a possible lag of 2-3 days of the VHE 
y-ray flare with respect to the lower frequencies. It should also 
be noted, that the X-ray flux is in general lower than during the 
MAGIC observations in 2006. 

In 2006-2007 there was also extensive whole earth blazar 
telescope (WEBT), very long baseline array (VLBA) and RXTE 
X-ray satellite monitoring of 3C 279 (Larionov et al. 2008). The 
VLBA data shows a component emerging from the core at MJD 
54063 ± 40 accompanied by the maxima in the 37 GHz light 
curve while the optical polarization data shows a rotation of the 
polarization angle ~ 300° starting around MJD 54115 with du- 
ration of ~ 2 months. There is coincident rotation of the polar- 
ization angle of the VLBI 43 GHz core and thus Larionov et al. 
(2008) conclude that the optical flares and rotation of the optical 
polarization angle take place in the VLBI knot as it moves down- 
stream of the core. The timing of the MAGIC detection (MJD 
54116) suggests that the emission of the VHE y-rays also hap- 
pens in this same emission region. This would be in agreement 
with previous results, that the y-ray flares observed by EGRET 
and Fermi would take place in the knots freshly emerging from 
the VLBI core (see Jorstad et al. 2001, Jorstad et al. 2010, and 
Agudo et al. 2011). This is sometimes denoted as 'far dissipa- 
tion' scenario, as the VLBI core is located tens of parsecs from 
the central engine (e.g. Marscher et al. 2010a). However, it has 
been argued (e.g. Tavecchio et al. 2010) that the size of the emis- 
sion region at this distance from the central engine is rather large 
which might be in contradiction with the day scale variability of 
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Fig. 5. Light curve in different energy bands during the 2007 ob- 
servation campaign. From top to bottom: VHE y-ray flux above 
150 GeV as measured by MAGIC, RXTE PCA flux in the 2-10 
keV range (Chatterjee et al. 2008), KVA R-band observations 
and REM infrared observations. For the MAGIC light curve the 
night by night flux is calculated assuming that 3C 279 always 
emits y-rays above 150 GeV. 



the MAGIC observation. This is further discussed in the follow- 
ing sections. 

In December 2008 - January 2009 the source was in a quies- 
cent state in the optical band. The polarization degree remained 
rather high (~30%) during the MAGIC observations and the po- 
larization angle was rotating smoothly from 60 to 90 degrees. 
The light curves are shown in Figure 6. In the X-ray band 
the source was at quiescent level, the flux being about 50% 
lower than in January 2007 (Abdo et al. 2010). In the y-ray 
band (Fermi) the source was active from mid-November 2008 
to the end of February 2009 and showed two major flares (first 
one peaking end of November (MJD 54795) and second mid 
February (MJD 54880)), but the core of the MAGIC observa- 
tions in 2009 January (MJD 54829-54856) took place in between 
the two flares, when the y-ray flux was at a relatively low level. 
Unfortunately, due to bad weather on the MAGIC site, no data 
could be taken in February, when Fermi observed the y-ray flare 
which was accompanied by a rotation of the optical polarization 
angle nor at the end of April when the source was flaring in the 
X-ray band (Abdo et al. 2010). 

The optical to VHE y-ray data of the three MAGIC observ- 
ing epochs is summarized in Fig. 7. In summary, although the 
detections at VHE y-ray energies correspond to high states in 
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Fig. 6. The optical R-band light curve (top), degree of polar- 
ization (middle) and polarization angle of 3C 279 from 2008 
December to 2009 February. The vertical lines show the time of 
MAGIC observations. 
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Fig. 7. SEDs of simultaneous optical, X-ray and y-ray data at 
the epochs of MAGIC observations. Blue: February 23, 2006, 
red: January 16, 2007, green January 21-February 1 2009. 



optical and X-rays, no multiwavelength correlations on short 
timescales seem to be present. 

5. Spectral Energy Distributions and Models 

In the following we discuss the spectral energy distributions 
(SEDs) of 3C 279 at three epochs, February 23, 2006, January 
16, 2007 and January 21-February 1, 2009. The first epoch cor- 
responds to the MAGIC discovery of the VHE y-ray emission 
from 3C 279 reported in Albert et al. (2008a). In addition to 



the multiwavelength data discussed there we show IR, optical 
and X-ray data simultaneous to the MAGIC observations from 
Bottcher et al. (2009). For January 16, 2007 and January 21- 
31, 2009 we show the quasi-simultaneous data described in the 
Section 3. 

The optical-IR data are important for SED modeling, since 
the steep continuum revealed by these measurements constrain 
the position of the synchrotron peak below the optical band in 
2006 and below the IR band in 2007. The X-ray spectral data 
are also quite constraining, since they allow us to define the low 
energy end of the inverse Compton component. 

The spectrum observed by MAGIC was corrected for inter- 
galactic absorption using the recent EBL model by Dominguez 
et al. (2010) adopted here throughout. In the original publication 
of the 2006 observations the EBL model by Primack et al. (2005) 
was used, which underestimates the IR background according to 
observations (e.g. Fazio et al. 2004; Chary et al. 2004). However, 
as discussed in Dominguez et al. (2010), the difference between 
the EBL models is smaller than the systematic uncertainties of 
the MAGIC data analysis. 

5.1. SED February 2006 

We try to reproduce the revised data with a simple one- 
zone leptonic emission model for FSRQs (details in Maraschi 
& Tavecchio 2003) considering the synchrotron and inverse 
Compton (IC) emission from a population of relativistic elec- 
trons in a single emission region (spherical, with radius R) in mo- 
tion with bulk Lorentz factor T at an angle 9 with the line of sight. 
The electron energy distribution is described by a smoothed 
broken-power law with normalization K (measuring the num- 
ber density of electrons at y = 1) extending from j\ to yj, with 
indices ri\ and m below and above the break at a yi,, respectively. 
The magnetic field with intensity B is supposed to be homoge- 
neous and tangled. The seed photons for the IC emission are both 
the synchrotron photons produced within the jet (SSC mecha- 
nism) and those outside the jet. We consider two cases, in which 
the IC process occurs within or outside the BLR. In the first case 
the high energy emission is dominated by comptonization of the 
UV photons of the BLR (EC/BLR). In the latter case the exter- 
nal radiation field is dominated by the IR thermal emission from 
the (putative) dusty torus (EC/IR, Ghisellini & Tavecchio 2008, 
Sikora et al. 2009). In the case of the BLR emission we assume 
that the clouds are located at a distance 7?blr from the central 
black hole. The resulting emission is modeled as a black body 
peaking at 10 15 Hz (Tavecchio & Ghisellini 2008) with lumi- 
nosity fixed to a fraction Tblr of the disk luminosity (fixed to 
L D = 2 x 10 45 erg s _1 , Pian et al. 1999 ). In the case of the torus 
emission, located at a distance we assume that a fraction 
tir = 0.5 of the disk luminosity is intercepted and re-radiated 
from dust as IR emission (again, we assume a black body spec- 
trum, with temperature T\g_ = 10 3 K, see Nenkova et al. 2008). 

In Figure 8 the blue line indicates the SED obtained applying 
the first model (EC/BLR), while the red line corresponds to the 
EC/IR case. The model parameters are listed in Table 4. 

The main difference between the two models (BLR and 
IR) is related to the different frequency of the target photons, 
determining the energy threshold above which the IC process 
is strongly suppressed due to the onset of the Klein-Nishina 
(hereafter KN) regime. Basically, this energy can be written as 
(Tavecchio & Ghisellini 2008, Ghisellini & Tavecchio 2009): 
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Fig. 8. The spectral energy distribution of 3C 279 on February 
23, 2006. The red symbols show the optical (filled circles), 
RXTE (bow-tie) Boettcher (2009) and MAGIC deabsorbed data 
(triangles) used to fit the two model curves: blue line assum- 
ing EC inside the BLR and red line EC outside the BLR. Also 
MAGIC observed data is shown (cyan). For comparison also 
historical data are shown: 1991 high state (gray: Hartman et al. 
1996) , 1993 low state (green: Maraschi et al. 1993) 1996 high 
state (orange, blue bow-tie: Ballo et al. 2002, Wehrle et al. 1998) 
and 2003 low state (magenta: Collmar et al. 2004). The dashed 
lines correspond to blackbody radiation from the IR torus (red) 
and BLR (blue). 



where 6 is the relativistic Doppler factor (6 — T for 6=1 /T) and 
v i5 is the frequency of the target photons (in units of 10 15 Hz). 
Therefore, if the IC process takes place within the BLR the emis- 
sion in the MAGIC band is strongly affected. The KN effect is 
more moderate if the IC scattering takes place outside the BLR. 
In this case the dominant photon frequency is v = 10 13 Hz cor- 
responding to Tjr = 10 3 K; hence £kn — 300 GeV. Accordingly, 
somewhat harder spectra are possible in the MAGIC band. 

The modeling of the VHE y-ray emission from 3C 279 
is particularly challenging (Albert et al. 2008a, Bottcher et al. 
2009). Indeed, due to the reduced IC scattering efficiency at high 
energy (e.g. Tavecchio & Ghisellini 2008) one zone emission 
models predict a low and rather soft emission in the VHE y-ray 
band. 

An additional difficulty for the BLR model is the possibility 
that VHE photons will be further reduced due to the internal 
yy — > e + e~ absorption with the UV photons of the BLR (Liu et 
al. 2008, Sitarek & Bednarek 2008, Tavecchio & Mazin 2008). 
In contrast, when the emission occurs outside the BLR, photons 
with E < 1 TeV are not substantially absorbed (above 1 TeV they 
will be absorbed interacting with the IR radiation of the torus). 
Summarizing, these observations are very difficult to reconcile 
with the one-zone framework. Both the EC/BLR and the EC/IR 
models can reproduce the data but require a rather large flux in 
the LAT band. Future detections of 3C 279 simultaneously in 
the LAT and VHE band will be crucial to confirm or rule out 
this possibility. 



5.2. SED January 2007 

At the time of the second MAGIC detection, in January 2007, 
3C 279 was in a brighter optical state but in a fainter X-ray state 
than at the discovery (February 23, 2006). Motivated by the dif- 
ficulties of the one-zone model discussed above (for 2007 the 
required GeV flux would be even higher than for 2006), we con- 
sider two-zone model and lepto-hadronic model for this dataset. 

5.2.1. Two-zone Model 

In"two zone" model, the emission from the optical up to the X- 
ray and y-ray bands derives from a different emission region than 
the VHE y-ray photons. As the VHE y-ray flare follows the op- 
tical one, we assume that the optical up to the X-ray and y-ray 
emission zone is closer to central engine than the VHE y-ray 
emission. Abdo et al. (2010) show that the y-ray emission re- 
gion is located on the parsec scale jet and we therefore model 
the spectral energy distribution assuming that the optical up to 
the X-ray and y-ray bands are emitted within the BLR while the 
VHE y-ray emission origins in a region outside the BLR. This 
also minimizes the effect of the KN regime and internal absorp- 
tion of y-rays. However, the modeling is valid also for emission 
regions further out in the jet as far as it is still located within IR 
torus to provide enough seed photons. 

Assuming two regions we are doubling the number of free 
parameters, therefore this scenario is far less constrained than 
one-zone models. Clearly the magnetic field and particle density 
will be lower and the size larger for the more external region. The 
parameters for the two regions are reported in Table 4 and the fit 
to 2007 data in Fig. 9. The two-zone model can better reproduce 
the MAGIC data, due to the possibility to shift the peak of the 
EC bump to higher energies. 

It is also interesting to speculate about the possible connec- 
tion between the energetic flare observed in the IR band and the 
VHE event, following the former by almost two days. In the 
framework of the two-zone model, an appealing picture assumes 
that the jet perturbation (possibly a shell) responsible for the IR 
flare (assumed to be produced inside the BLR), travels and even- 
tually reaches the region where VHE photons can be more eas- 
ily produced and can escape the emission region. From the delay 
between the two events we can infer the distance between the lo- 
cation of the infrared to optical synchrotron and VHE emission 
zones, d - ct\^T 2 ~ 2 x 10 18 cm, where we assume that the 
perturbations travels with a Lorentz factor F ^ 20, the value as- 
sumed in the radiative model. Interestingly, this distance is very 
well consistent with the assumed size of the VHE y-ray emission 
zone. 

5.2.2. Lepto-hadronic model 

Although hadronic models in general are not favored for lu- 
minous quasars (Sikora et al. 2009), they might be viable for 
3C 279 as it does not have as hard X-ray spectrum (F < 
0.5, for 3C 279 F=0.66 Ballo et al. 2002) as other luminous 
quasars. Furthermore, models taking into account both leptonic 
and hadronic processes could be viable for describing the spec- 
tral energy distribution of quasars. Therefore, we also fit the data 
with a lepto-hadronic model. The model will be described in full 
detail in a forthcoming paper (Riiger et al. in prep.). It com- 
prises a non-thermal proton and electron distribution injected 
into the radiation volume as a power law with lower and up- 
per bound (y mm and ymax respectively). The particles may ra- 
diate through a number of channels: The electrons may radiate 
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Fig. 9. The spectral energy distribution of 3C 279 on January 
2007 16th (MJD 541 16). The red symbols show the REM (MJD 
54114.4, triangles), KVA (filled circle), RXTE (bow-tie) and 
MAGIC (deabsorbed, triangles) data used to fit the two-zone 
model (see text). Blue line shows the emission from the zone 
inside BLR and red line the emission from outside the BLR. 
Additionally REM data from MJD 54113.3 (Maximum of the 
REM lightcurve) cyan and MJD 541 17.3 (black) are shown. The 
dashed lines correspond to blackbody radiation from the IR torus 
(red) and BLR (blue). The historical data as in Fig. 8. 



through the above described mechanisms of synchrotron emis- 
sion and inverse Compton scattering off synchrotron photons. 
The protons will also emit synchrotron radiation. Additionally 
they take part in photo-hadronic processes. Those are modeled 
by the cross-sections described in Kelner & Aharonian (2008). 
In the photo-hadronic processes tt° and 1& are produced. The 
7T° decay into y-rays may eventually start a pair cascade in the 
source if y-rays are emitted in the opticallythick regime. The 
7r ± decays into electrons and positrons which again emit syn- 
chrotron radiation. The time evolution is solved self-consistently 
with a scheme described in Riiger et al. (2010) for the pure lep- 
tonic case. 

In this model the number of free parameters is comparable 
to one-zone SSC models. We introduce four extra parameters to 
describe the proton population. Since we assume identical mini- 
mum Lorentz factor and spectral index for electrons and protons, 
we end up with a total of nine free parameters. Unlike the one- 
zone SSC model the variability patterns are far more compli- 
cated. In this model the lower limit of variability is still given by 
the beamed light crossing time, which is of the order of 4 days. 
Since electrons and protons may burst differently and at different 
times, the determination of the variability pattern is complicated. 
If electron and proton density are increased simultaneously, the 
electrons produce the well-known soft lag pattern in the electron 
synchrotron radiation and the inverse Compton radiation. In the 
hadronic component this is also true for the synchrotron radia- 
tion. For the n production intermediate secondary energies peak 
earlier. 

To model the data (fit shown in Fig. 10) we used electron and 
proton spectra with a y mm = 150 and a spectral index s = 2.2. 
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Fig. 10. The observed spectral energy distribution of January 
2007 modeled with the lepto-hadronic model. The overall fit is 
shown with red line with the following components at high en- 
ergies: the synchrotron radiation of pair creation electrons and 
positrons cascaded down from the optically thick regime (ma- 
genta dotted line), synchrotron radiation from positrons from the 
pion decay (green dashed line), synchrotron radiation from elec- 
trons from the pion decay (blue dashed line), inverse Compton 
scattering (cyan dot dashed line) and proton synchrotron emis- 
sion (black double dashed line). The VHE y-ray emission is 
mostly sum of the three first components while in X-rays the 
main contribution comes from the inverse Compton scattering 
(like in the purely leptonic models). The low energy bump is 
produced by the electron synchrotron radiation. 



For the protons we find ymax = 10 9 , while the electrons have a 
much lower ymax = 5 x 10 4 . Spectral breaks are calculated self- 
consistently from the loss processes within the radiation zone. 
The magnetic field is assumed to be 0.025 G with a source radius 
of 5.3 x 10 17 cm. The doppler factor is 42. The energy density of 
the electrons is one order of magnitude higher than the magnetic 
field energy density, while the proton energy density is already 
six orders of magnitude higher. Even though this is not in equi- 
librium this might still be confined, as the proton gyroradius at 
highest energies is of the order of the blob radius. 

The parameters suggest that the radiation zone may be close 
to the end of magnetic confinement zone, ~1000 Schwarzschild 
radii from the supermassive black hole. However, the model it- 
self is not limited to that region, but may be used on every point 
along the jet axis. Although internal pair absorption is included, 
external photons producing pairs are not included in the model 
and therefore if the emission region was inside the BLR the ob- 
served flux would be reduced. 



5.3. SED January 2009 

For the 2009 observations we present a SED built using multi- 
wavelength data nearly simultaneous to the MAGIC data taken 
in the period 21-31 January, which yielded in an upper limits 
in the VHE y-ray band, y-ray data were derived from Fermi 
LAT and averaged over the same period (Fermi data were not 
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Fig. 11. The spectral energy distribution of 3C 279 on January 
2009: MAGIC (taken in the period 21-31 January, red arrows: 
deabsorbed, cyan arrows: observed), LAT (averaged over the 
same period, red pentagons), XRT and UVOT (February 1, 2009 
red squares) and KVA (January 25, 2009 red filled circle). The 
SED is modeled assuming the emission region inside the BLR. 
The historical data as in Fig. 8. The dashed line correspond to 
blackbody radiation from the IR torus (red). 

available at the previous epochs) and XRT and UVOT data of 
February 1, 2009 were used. The source was in a rather low state 
at all available bands. The SED can be modeled quite satisfacto- 
rily with a standard one zone model assuming the emission re- 
gion inside the BLR (Fig. 1 1) with typical parameters (reported 
in Table 4). 

6. Summary and Conclusions 

In this paper we report the MAGIC observations of 3C 279 
in January 2007 and from December 2008 to April 2009. The 
January 2007 observations yielded a detection on the night of 
January 16, when the source was in a high optical state. Like in 
February 2006 when the source was discovered, in January 2007 
the source was detected on one night only. However, in January 
2007 MAGIC observed the source also in the previous and the 
following night and these observations resulted in non-detection, 
which confirms the day-scale variability of the source. 

In February 2006 and January 2007 the source was in a high 
optical-IR state and was also showing a rather high state in X- 
rays. However, no major increase in the optical or X-ray flux was 
seen simultaneously to the VHE y-ray flare. Unlike in February 
2006 and January 2007, in January 2009 the source was observed 
by MAGIC during a low optical-X-ray state. The existence of 
LAT data for this epoch shows that the source was relatively 
faint also at GeV energies. Consistently, the MAGIC observa- 
tions resulted in upper limits. 

In the spring of 2009 Fermi detected a fast y-ray flare from 
3C 279 which was coincident with a rotation of the optical po- 
larization angle of 180 degrees that had duration of ~ 20 days. 
In early 2007 such rotation of the polarization angle also took 
place (Larionov et al. 2008), but with rotation of ~ 360 degrees, 



the duration of rotation ~2 times longer and the direction of the 
rotation opposite to what was seen in 2009. The MAGIC detec- 
tion in January 2007 took place in the beginning of the rotation 
of optical polarization angle, thus confirming that rotation of op- 
tical polarization angle in 3C 279 are recurrently accompanied 
with y-ray flares. This seems to be common in FSRQs, the same 
behavior has been observed in PKS 15 10-089 in 2009 (Marscher 
et al. 2010b). 

However, the patterns in 3C 279 in January 2007 and 
February 2009 are clearly different, which might suggest that 
different mechanisms for producing the rotation of the polar- 
ization angle could be in action (e.g. for 2009 the rotation 
could be produced by turbulence, see e.g. D'Arcangelo et al. 
2007). Therefore it is not possible to draw definite conclusions 
on the connection between the rotation of polarization angle with 
y-ray flares either. 

We find that the spectral energy distribution of 3C 279 in 
2007 January cannot be explained by the standard one-zone 
SSC+EC model, where the emission region is located inside the 
BLR which provides the seed photons for inverse Compton scat- 
tering. Instead we explore a two-zone model, where the VHE 
y-ray emitting region is located just outside the BLR, while the 
standard optical-to-X-ray and y-ray emitting region is still in- 
side the BLR region and a lepto-hadronic model, which both fit 
the data reasonably well. Also other possible scenarios for the 
emission exists. Sitarek & Bednarek (2010) suggested a leptonic 
cascade model. Several authors have also suggested that blazar 
emission zone would be located at parsec scale distances from 
the black hole: there y-rays could be produced via up-scattering 
of infrared photons from the hot dust (Blazejowski et al. 2000, 
Sikora et al. 2009). It has also been suggested that the relativistic 
jet could drag part of the BLR clouds to large distances from the 
central engine (at least in radio galaxies: Arshakian et al. 2010, 
Leon-Tavares et al. 2010a) and in principle these photons could 
serve as external seed photons for inverse Compton scattering 
(Leon-Tavares et al. 201 1). The multi wavelength behavior of the 
source in January 2007 can also be interpreted in favor of emis- 
sion zone beyond VLBI core and in such case multi emission 
shell model suggested by Marscher et al. (2010) could be viable. 

After the first detection of 3C 279 by MAGIC, two other 
flat spectrum radio quasars have been detected in the VHE y- 
ray band, PKS 1510-089 (z = 0.36, Wagner et al. 2010) and 
PKS 1222+21 (z = 0.432, Mariotti et al. 2010, Aleksic et al. 
2011). In all the cases the emission of the VHE y-ray radiation 
poses problems to the standard models, as discussed in this paper 
and in Aleksic et al. (2011). Future multiwavelength observa- 
tions, including higher sensitivity observations with MAGIC II, 
will be needed to improve our understanding of the high-energy 
emission from FSRQs. 
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